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Abstract 

Background: China is a country that is most seriously affected by hemorrhagic fever with renal syndrome (HFRS) 
with 90% of HFRS cases reported globally. At present, HFRS is getting worse with increasing cases and natural foci 
in China. Therefore, there is an urgent need for monitoring and predicting HFRS incidence to make the control of 
HFRS more effective. In this study, we applied a stochastic autoregressive integrated moving average (ARIMA) 
model with the objective of monitoring and short-term forecasting HFRS incidence in China. 

Methods: Chinese HFRS data from 1975 to 2008 were used to fit ARIMA model. Akaike Information Criterion (AIC) 
and Ljung-Box test were used to evaluate the constructed models. Subsequently, the fitted ARIMA model was 
applied to obtain the fitted HFRS incidence from 1978 to 2008 and contrast with corresponding observed values. 
To assess the validity of the proposed model, the mean absolute percentage error (MAPE) between the observed 
and fitted HFRS incidence (1978-2008) was calculated. Finally, the fitted ARIMA model was used to forecast the 
incidence of HFRS of the years 2009 to 201 1. All analyses were performed using SAS9.1 with a significant level of p 
< 0.05. 

Results: The goodness-of-fit test of the optimum ARIMA (0,3,1) model showed non-significant autocorrelations in 
the residuals of the model (Ljung-Box Q statistic = 5.95,P = 0.31 13). The fitted values made by ARIMA (0,3,1) model 
for years 1978-2008 closely followed the observed values for the same years, with a mean absolute percentage 
error (MAPE) of 12.20%. The forecast values from 2009 to 201 1 were 0.69, 0.86, and 1.21 per 100,000 population, 
respectively. 

Conclusion: ARIMA models applied to historical HFRS incidence data are an important tool for HFRS surveillance in 
China. This study shows that accurate forecasting of the HFRS incidence is possible using an ARIMA model. If 
predicted values from this study are accurate, China can expect a rise in HFRS incidence. 



Background 

Hemorrhagic fever with renal syndrome (HFRS), or epi- 
demic hemorrhagic fever (EHF) is an acute viral syn- 
drome caused by infection with one of hantaviruses. 
HFRS is an important infectious disease in developing 
countries. In China, HFRS is caused mainly by 2 types 
of hantaviruses, Hantaan virus (HTNV) and Seoul virus 
(SEOV), each of which has coevolved with a distinct 
rodent host. HTNV is associated with Apodemus agrar- 
ius, whereas SEOV, which causes a less severe form of 
HFRS, is associated with Rattus norvegicus [1]. In 
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hantavirus -endemic areas, HFRS is most common 
among farmers and others who may have close contact 
with excreta of infected rodents [2,3]. In mainland 
China, HFRS remains a serious public health problem 
with approximately 20,000-50,000 human cases reported 
annually, approximately 90% of the total cases world- 
wide [4-6]. Currently, HFRS is endemic in 28 of 31 pro- 
vinces in mainland China [4,7]. 

In response to the spread of HFRS in China, the Chi- 
nese Center for Disease Control and Prevention 
designed a surveillance system for HFRS and created 
educational programs for the general public. However, 
the impact of control efforts remains difficult to mea- 
sure due to the inherent complexities of HFRS as a dis- 
ease: multiple viral strains with identified genetic 
polymorphisms, complex disease manifestation, diverse 
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animal reservoirs, and multiple routes of transmission [8]. 
Infectious diseases have certain characteristic features that 
lead themselves to modeling, such as: speed of pathogen 
variation, accumulation of susceptible hosts, and environ- 
mental indices [9]. Thus, epidemic modeling and forecast- 
ing can be essential tools to prevent and control HFRS. 
Recently, statistical methods including linear regression 
[10-12], correlation coefficients [13], grey swing model 
[14], back propagation artificial neural network model [15] 
have been used for prediction of HFRS incidence. The var- 
iation of HFRS incidence, which is influenced and con- 
strained by diversified factors, is characterized by tendency 
and randomicity. These statistical tools are inappropriate 
for analyzing the randomicity of HFRS. Autoregressive 
integrated moving average (ARIMA) models, which take 
into account changing trends, periodic changes, and ran- 
dom disturbances in time series, are very useful in model- 
ing the temporal dependence structure of a time series. In 
epidemiology, ARIMA models have been successfully 
applied to predict the incidence of infectious diseases, 
such as influenza mortality [16], malaria incidence [17], as 
well as other infectious diseases [18,19]. This study aimed 
to develop a univariate time series model for the HFRS 
incidence; specifically, a stochastic ARIMA model, for 
short-term forecasting of HFRS incidence (per 100,000 
population) in China. 

Methods 

Materials 

Chinese HFRS incidence data from 1975 to 2008 was 
obtained from the Chinese Center for Disease Control 
and Prevention. All HFRS cases were initially diagnosed 
by clinical symptoms. Patient blood samples were also 
collected and sent to local Centers for Disease Control 
and Prevention (CDC) laboratories for serological con- 
firmation. Finally, data were collected by case number 
according to the sampling results. There might be 
admission rate bias in the disease report, but this has 
been reduced as much as possible. In China, HFRS is a 
nationally notifiable disease and hospital physicians 
must report every case of HFRS to the local health 
authority within 12 hours. Local health authorities later 
report monthly HFRS case totals to higher the national 
level CDC for surveillance purposes. Due to mandatory 
reporting, it is believed that the degree of compliance in 
disease notification over the study period was consistent. 

We used the Box-Jenkins approach to ARIMA (p, d, 
q) modeling of time series [20]. This model-building 
process is designed to take advantage of associations in 
the sequentially lagged relationships that usually exist in 
periodically collected data [21]. The following were the 
parameters selected when fitting the ARIMA model: p, 
the order of autoregression; d, the degree of difference; 
q, the order of moving average. 



The annual data used in this study did not show sea- 
sonal pattern, so the series was differenced at the non- 
seasonal level to induce stationarity. Autocorrelation 
function (ACF) graph and Partial autocorrelation func- 
tion (PACF) graph were used to identify the order of 
moving average (MA) and autoregressive (AR) terms 
included in the ARIMA model. Estimates of the model's 
parameters were obtained by the conditional least 
squares method. Diagnostic checking including residual 
analysis and the Akaike Information Criterion (AIC) was 
used to compare the goodness-of-fit among ARIMA 
models. The Ljung-Box test was used to measure the 
ACF of the residuals. In addition, we used the mean 
absolute percentage error (MAPE) and fitting effect dia- 
gram to assess forecast accuracy. 

MAPE = -J2\( x t~ xt)/xt\ , where x t and x t denote 

n t=i 

observed and fitted values at time point t. The MAPE 
value was calculated based on observed values and fitted 
values from 1978 to 2008. A lower MAPE value indi- 
cates a better fit of the data. Finally, the fitted ARIMA 
model was used for short-term forecasting of HFRS inci- 
dence for years 2009 to 2011. All analyses were per- 
formed using SAS9.1 with a significant level of p < 0.05. 

Ethical review 

The present study was reviewed by the research institu- 
tional review board of Shandong University and the 
China CDC, and found that utilization of disease sur- 
veillance and meteorological data did not require over- 
sight by an ethics committee. 

Results 

Temporal analysis 

From 1975 to 1986, the HFRS incidence in China rose 
regularly with a peak in 1986 of 11.06 cases per one 
hundred thousand population. After 1986, the incidence 
descended sharply with a dramatic fluctuation until 
2008 (Figure 1). The lowest incidence could be seen in 
2008, 0.68 per one hundred thousand. 

Model identification 

According Figure 1, the series showed a non-stationary 
mean, so we stabilized the mean of HFRS incidence by 
taking both second and third order differences. All 
further statistical procedures were performed on the 
transformed HFRS incidence. Based on the distribution 
characteristics, we conducted five models, ARIMA(0, 2, 
1), ARIMA(1,2,1), ARIMA(0,3,1), ARIMA(1, 3, 1), and 
ARIMA(2, 3, 1). Of all the models tested, the ARIMA 
(0,3,1) model was the best fit for the data (Table 1). The 
transformation series by taking third-order differences is 
shown in Figure 2. The plots of ACF and PACF (Figure 
3) described the temporal dependence structure in 
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Figure 1 The incidence of HFRS in China from 1975 to 2008 



HFRS incidence. The slow decay in the PACF, associated 
with a ACF cutoff at lag 1 suggested a MA(q = 1). 

Model diagnosis 

The parameter estimates for the optimum ARIMA(0,3,1) 
model are shown in Table 2. The model's fitted (1975- 
2008)and predicted values (2009-2011) are presented in 
Figure 4. The MAPE value was 12.20%. The forecast 
values of the years 2009, 2010, and 2011 were 0.69, 0.86, 
and 1.21 per 100,000 population, respectively. 

Discussion 

Time series analysis of surveillance data on incidence of 
various infections is very helpful in developing hypoth- 
eses to explain and anticipate the dynamics of the 
observed phenomena and subsequently in the establish- 
ment of a quality control system and reallocation of 
resources [22]. ARIMA model is one of the most widely 
used time-series forecasting techniques because of its 
structured modeling basis and acceptable forecasting 
performance [23]. In this paper, we applied an ARIMA 
(p, d, q) model to analyze the surveillance data of HFRS 
in China. Disease monitoring by public health 



Table 1 Comparisons of tested models 


Model 


Ljung-Box Q statistic 


P value 


AIC 


ARIMA(0,2,1) 


5.78 


0.3869 


117.7200 


ARIMA(1,2,1) 


7.97 


0.0488 


119.6435 


ARIMA(0,3,D 


5.95 


0.3113 


116.7075 


ARIMA(1,3,D 


8.46 


0.0453 


118.6647 


ARIMA(2,3,D 


9.33 


0.0260 


120.1425 



department entails ongoing data collecting, processing, 
and updating. However, the national level China CDC is 
the appropriate level of organization for the implemen- 
tation of an ARIMA predictive model, because reported 
data is continually received and updated. We found that 
model predictions are further improved by the assured 
availability of the Health Department data. In this study, 
we have obtained an ARIMA model that closely fits 
HFRS incidence in China. The autoregression and mov- 
ing average parameters of our model imply the inci- 
dence of HFRS in a month can be estimated by the 
residual occurring one month prior. According to the 
results above, the conducted model is reliable with a 
high validity. Once a satisfactory model has been 
obtained, it can be used to forecast expected numbers of 
cases for a given number of future time intervals [24]. 
Thus, the fitted ARIMA(0,3,1) model can be used to 
predict the next three years' HFRS incidence in China. 
The forecast results suggest that the HFRS incidence in 
China will experience a slight growth in the next three 
years (2009-2011). A rise in the number of HFRS inci- 
dence may also result from an increase in the number 
and size of natural foci [25], climate change, especially 
the increase of mean temperature [26,27]. Therefore, 
knowledge of HFRS forecasts is necessary to prompt 
health departments to strengthen surveillance systems 
and reallocate resources in anticipation of increasing 
HFRS incidence. 

Several studies have used ARIMA model to fit and 
predict changing trends in infectious disease. Luz et al 
applied an ARIMA(2,0,0)x(l,0,0) 12 model to predict 
dengue incidence in Rio de Janeiro [18] and found that 
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Figure 2 Third-order differences of HFRS incidence(per 100,000 pop.) 

V " J 

ARIMA models were useful tools for monitoring dengue applied an ARIMA model to monthly incidence of 
incidence. Earnest et al indicated that ARIMA models HFRS in Linyi City, China to predict HFRS incidence, 
provided useful tools for administrators and clinicians in and found that the ARIMA model could be used to pre- 
planning for real-time bed capacity during infectious diet HFRS incidence with high predictive precision in 
diseases outbreaks such as SARS [28]. Li et al have the short-term [29]. In the present study, we further 
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Figure 3 Autocorrelation(acf / top) and partial autocorrelation(pacf, bottom) functions of third-order differences HFRS incidence. Dotted 
line: 95% confidence intervals. 
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Table 2 Parameter for the final ARIMA(0,3,1 ) model 



Parameter 


Coefficient 


Standard error 


t statistic 


P value 


MA1 


0.9675 


0.0599 


16.14 


< 0.0001 



confirmed the consensus that ARIMA model is a useful 
tool in monitoring and predicting changing trends in 
infectious diseases. 

To the best of our knowledge, this is the first study to 
apply ARIMA model to fit the HFRS incidence in China 
with as many as 34 observations at year level. Some pre- 
vious studies [30,31] in China also used ARIMA model 
to fit and forecast HFRS incidence of some regions, but 
they had the same problem that the number of observa- 
tions was not enough, which led to the instability of 
their forecast results. In order to conduct a stable and 
effective ARIMA model, we have to collect at least 30 
observations [32]. Thus, parameter estimates of the 
fitted model would be more robust. The longer the ser- 
ies, the better; however, the series should not extend so 
far into the past as to include periods during which a 
different case definition was applied or in which any 
other reporting artifact resulted in a mean number of 
cases per interval that differs from the mean of recent 
intervals. As mentioned above, for adequate ARIMA 
modeling, a time series should be stationary with respect 
to mean and variance. If the mean increases or 
decreases over time, or if the variance does, the series 
may need to be transformed to make it stationary, 
before being modeled. Otherwise, the prediction effect 
of the model will be poor. 



In order to improve the model, updating the forecasts 
is very important. A model without seasonal terms will 
need to be updated frequently. Confidence intervals that 
widen rapidly as time increase from the starting point of 
the forecasts also indicate a model that needs frequent 
updating. Generally speaking, there are two ways to 
implement the updating. The model can be reapplied to 
the original series with extra observations added at the 
end to give forecasts based on a later starting point. 
Alternatively, a new model can be fitted to the longer 
series. This is probably preferable, since fitting a model 
is quick, especially when the old model is used as a 
guide, and it makes better use of the additional 
observations. 

Some limitations of this study also need to be taken 
into account when interpreting the results. In this study, 
the interval of HFRS incidence is one year, so we could 
not analyze its seasonal characteristic. In further study, 
we would use monthly data to predict HFRS incidence 
in order to get seasonal pattern and higher predictive 
precision. In addition, the data are from a passive sur- 
veillance system, the possible biases in disease reporting 
and potential underreporting of HFRS cases might influ- 
ence the precision of our analysis. 

Conclusion 

There is an urgent need for monitoring and predicting 
HFRS incidence to reduce the substantial morbidity and 
mortality caused by this disease [33]. ARIMA models 
applied to historical HFRS incidence data are an impor- 
tant tool for HFRS surveillance. Accurate forecasting of 
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Figure 4 Black dots: observed values of HFRS incidence, for the period 1975-2008. Red solid line: ARIMA(0,3J) model's fitted curve(1978- 
2008). Red dots: ARIMA(0,3,D model's fitted values(1 978-2008). Dashed lines: 95% confidence intervals. 
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the incidence of HFRS is possible. Our modeling 
approach can be used to monitor and predict HFRS 
incidence in China. The ARIMA model could be used 
to optimize HFRS prevention by providing estimates on 
HFRS incidence trends in China. 

Acknowledgements 

This study was supported by the social welfare research special program of 
the Ministry of Science and Technology, China (2003DIA6N009) and Special 
Infectious Diseases Program of the Ministry of Science & Technology, China 
(Grant No.2008ZX1 0004-010). 

Author details 

1 National Institute for Communicable Disease Control and Prevention, China 
CDC, Beijing, 102206, PR China. 2 State Key Laboratory for Infectious Diseases 
Prevention and Control, Beijing, 102206, PR China, department of 
Epidemiology and Health Statistics, School of Public Health, Shandong 
University, Jinan, 250012, Shandong Province, PR China. 4 Chinese Center for 
Disease Control and Prevention, Beijing, 102206, PR China. 

Authors' contributions 

QL, XL, BJ and WY conceived the study, undertook statistical analysis and 
drafted the manuscript. XL and BJ assisted with data collection and 
statistical analysis. All authors contributed to the writing of the manuscript 
and approved the submitted version of the manuscript. 

Competing interests 

The authors declare that they have no competing interests. 

Received: 18 November 2010 Accepted: 15 August 2011 
Published: 15 August 2011 

References 

1. Fang LQ, Zhao WJ, Vlas SJ, Zhang WY, Liang S, Looman CWN, Yan L, 
Wang LP, Ma JQ, Feng D, Yang H, Cao WC: Spatiotemporal dynamics of 
hemorrhagic fever with renal syndrome, Beijing, People's Republic of 
China. Emerging Infectious Diseases 2009, 15:2043-2045. 

2. Vapalahti K, Paunio M, Brummer-Korvenkontio M, Vaheri A, Vapaahti O: 
Puumala virus infections in Finland:increased occupational risk for 
farmers. Am J Epidemiol 1999, 149(12):1 142-1 151. 

3. Glass GE, Childs JE, Korch GW, LeDuc JW: Association of intraspecific 
wounding with hantaviral infection in wild rats(Rattus norvegicus). 
Epidemiol Infect 1988, 101:459-472. 

4. Zhang YZ, Xiao DL, Wang Y, Wang HX, Sun L, Tao XX, Qu YG: The 
epidemic characteristics and preventive measures of hemorrhagic fever 
with renal syndrome in China. Chin J Epidemiol 2004, 25(6):466-469. 

5. Ulrich R, Hjelle B, Pitra C, Kuger DH: Emerging viruses: the case 
'hantavirus'. Intervirology 2002, 45(4):31 8-322. 

6. Fang LQ, Yan L, Liang S, Vlas SJ, Feng D, Han X, Zhao W, Xu B, Bian L, 
Yang H, Gong P, Richardus JH, Cao WC: Spatial analysis of hemorrhagic 
fever with renal syndrome in China. BMC Infect Dis 2006, 6:77. 

7. Yan L, Fang LQ, Huang HG, Zhang LQ, Feng D, Zhao WJ, Zhang WY, Li XW, 
Cao WC: Landscape elements and Hantaan virus-related hemorrhagic 
fever with renal syndrome, People's Republic of China. Emerg Infect Dis 
2007, 13(9):1 301 -1306. 

8. Zhang Y: The epidemiological research status and problems and 
prospects of hemorrhagic fever with renal syndrome in China. Chin J 
Vector Bio & Control 2002, 13(2):85-88. 

9. Guan P, Huang DS, Zhou BS: Forecasting model for the incidence of 
hepatitis A based on artificial neural network. World J Gastroenterol 2004, 
10(24)3579-3582. 

10. Wang YJ, Zhao TQ, Wang P, Li SQ, Huang Z, Yang GQ, Li XY, Liu B: 
Applying linear regression statistical method to predict the epidemic of 
hemorrhagic fever with renal syndrome. Chin J Vector Bio & Control 2006, 

17(4)333-334. 



11. Olsson GE, Hjertqvist M, Lundkvist A, Hornfeldt B: Predicting high risk for 
human hantavirus infections, Sweden. Emerg Infect Dis 2009, 
15(1):104-106. 

1 2. Bi P, Wu XK, Zhang FZ, Parton KA, Tong SL: Seasonal rainfall variability, the 
incidence of hemorrhagic fever with renal syndrome, and prediction of the 
disease in low-lying areas of China. Am J Epidemiol 1998, 148(3):276-281. 

13. Clement J, Vercauteren J, Verstraeten WW, Ducoffre G, Barrios JM, 
Vandamme AM, Maes P, Ranst MV: Relating increasing hantavirus 
incidences to the changing climate: the mast connection. International J 
of Health Geographies 2009, 8:1. 

14. Guo LC, Wu W, Guo JQ, Wang P, Zhou BS: Appling grey swing model to 
predict the incidence trend of hemorrhagic fever with renal syndrome 
in Shenyang. Journal of China Medical University 2008, 37(6):839-842. 

15. Wu ZM, Wu W, Wang P, Zhou BS: Prediction for incidence of hemorrhagic 
fever with renal syndrome with back propagation artificial neural 
network model. Chin J Vector Bio & Control 2006, 1 7(3):223-226. 

16. Reichert TA, Simonsen L, Sharma A, Pardo SA, Fedson DS, Miller MA: 
Influenza and the winter increase in mortality in the United States, 
1959-1999. Am J Epidemiol 2004, 160(5):492-502. 

17. Gaudart J, Toure O, Dessay N, Dicko AL, Ranque S, Forest L, Demongeot J, 
Doumbo OK: Modelling malaria incidence with environmental 
dependency in a locality of Sudanese savannah area, Mali. Malaria 
Journal 2009, 8:61. 

18. Luz PM, Mendes BV, Codeco CT, Struchiner CJ, Galvani AP: Time series 
analysis of dengue incidence in Rio de Janeiro, Brazil. Am J Trop Med Hyg 
2008, 79(6):933-939. 

19. Yi J, Du CT, Wang RH, Liu L: Applications of multiple seasonal 
autoregressive integrated moving average(ARIMA) model on predictive 
incidence of tuberculosis. Chinese Journal of Preventive Medicine 2007, 
41(2):118-121. 

20. Box GEP, Jenkins GM: Time series analysis: forecasting and control. San 
Francisco: Holden Day; 1976, 181-218. 

21. Akhtar S, Rozi S: An autoregressive integrated moving average model for 
short-term prediction of hepatitis C virus seropositivity among male 
volunteer blood donors in Karachi, Pakistan. World J Gastroenterol 2009, 
15(1 3):1 607-1 61 2. 

22. Kuhn L, Davidson LL, Durkin MS: Use of poisson regression and time 
series analysis for detecting changes over time in rates of child injury 
following a prevention program. Am J Epidemiol 1994, 140(1 0):943-955. 

23. Wong J, Chan A, Chiang YH: Time series forecasts of the construction 
labour market in Hong Kong: the Box-Jenkins approach. Construction 
Management and Economics 2005, 23(9):979-991. 

24. Allard R: Use of time-series analysis in infectious disease surveillance. 
Bulletin of the World Health Organization 1 998, 76(4):327-333. 

25. Wang XF, Wang MW, Sun H: Epidemiological analysis of hemorrhagic 
fever with renal syndrome in China from 2004 to 2005. Disease 
Surveillance 2007, 22(5):307-309. 

26. Bi P, Tong SL, Donald K, Parton K, Ni JF: Climatic, reservoir and 
occupational variables and the transmission of hemorrhagic fever with 
renal syndrome in China. Int J of Epidemiol 2002, 31:189-193. 

27. Clement J, Vercauteren J, Verstraeten WW, Ducoffre G, Barrios JM, 
Vandamme AM, Maes P, Ranst MV: Relating increasing hantavirus 
incidences to the changing climate: the mast connection. International 
Journal of Health Geographies 2009, 8:1. 

28. Earnest A, Chen Ml, Ng D, Leo YS: Using autoregressive integrated 
moving average(ARIMA) models to predict and monitor the number of 
beds occupied during a SARS outbreak in a tertiary hospital in 
Singapore. BMC Health Services Research 2005, 5:36. 

29. Li XJ, Kang DM, Cao J, Wang JZ: A time series model in incidence 
forecasting of hemorrhagic fever with renal syndrome. Journal of 
Shandong University (Health Sciences) 2008, 46(5):547-549. 

30. Wu W, Guan P, Guo JQ, Zhou BS: Comparison of GM(1,1) gray model and 
ARIMA model in forecasting the incidence of hemorrhagic fever with 
renal syndrome. Journal of China Medical University 2008, 37(1 ):52-55. 

31. Chen Y, Bai S, Chen HZ, Sun BJ, Wei WJ, Huang M, Wang P: Fitting 
research on ARMA model in the prediction of incidence trend of 
hemorrhagic fever with renal syndrome. Modern Preventive Medicine 2008, 
35(8):1414-1415. 



Liu et al. BMC Infectious Diseases 2011, 11:218 
http://www.biomedcentral.eom/1 471 -2334/1 1/218 



Page 7 of 7 



32. Gao HX compilation: SAS System-SAS/ETS Software Manual. Beijing: China 
Statistics Press;83. 

33. Chen HX, Luo CW: Surveillance of hemorrhagic fever with renal 
syndrome in China. Chin J Epedimiol 2002, 23(1):63-66. 

Pre-publication history 

The pre-publication history for this paper can be accessed here: 
http://www.biomedcentral.eom/1471-2334/1 1/218/prepub 



doi:1 0.1 1 86/1 471-2334-1 1-218 

Cite this article as: Liu et al:. Forecasting incidence of hemorrhagic 
fever with renal syndrome in China using ARIMA model. BMC Infectious 
Diseases 2011 11:218. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 

Submit your manuscript at (^\ RioMM i r pnt ^\ 

www.biomedcentral.com/submit Blomea ^ enirai 



